Comment on 'Distorted perovskite with configuration as a frustrated spin system' 
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Interesting magnetic structure found experimentally in a series of manganites RMnOa, where 
R is a rare earth, was explained, ostensibly, by Kimura et al, in terms of a classical Heisenberg 
model with an assumed set of exchange parameters, nearest-neighbor and next-nearest-neighbor 
interactions, J„. They calculated a phase diagram as a function of these parameters for the ground 
state, in which the important "up-up-down-down" or E-phase occurs for finite ratios of the J„. In 
this Comment we show that this state does not occur for such ratios; the error is traced to an 
incorrect method for minimizing the energy. The correct phase diagram is given. We also point out 
that the finite temperature phase diagram presented there is incorrect. 

PACS numbers: 75.10.Hk,75.30.Kz,75.47.Lx 



Interesting magnetic phases that had been observed in 
a series of manganites RMnOa , with R=holmium in par- 
ticular, were explained in the cited paper [l] on the basis 
of a Heisenberg Hamiltonian with classical spins on a 2-d 
rectangular lattice with ferromagnetic nearest-neighbor 
exchange parameter, Ji, and two different antiferromag- 
netic next-nearest- neighbor parameters, J2 and J3. (The 
difference is due to the asymmetric positions of the oxy- 
gens resulting from Jahn- Teller distortion). By consid- 
ering a small set of states including the unusual E-type, 
or "up-up-down-down" state, they obtained a phase dia- 
gram for the ground state, and at finite temperature T. 
We show that the ground state phase diagram is incor- 
rect, and point out that the finite-T case is also incor- 
rect, the correct results differing qualitatively from those 
presented there. The importance of these conclusions is 
heightened by the fact that detailed results of the model 
calculation are accepted in 

Circa I960 it was shown that guessing various spin 
states ("reasonable", or experimentally determined) as 
ground states of the classical Heisenberg energy. 



H — ^ ' Jij Sj 



Sj , where 



Sf = 1, 



(1) 



with some assumed set of exchange parameters Jij , and 
comparing their energies to arrive at a phase diagram, is 
fraught with danger, often giving incorrect results. (The 
spin lengths have been absorbed into the J's.) Further- 
more, for spins on a Br avals lattice, it was shown Q that 
the ground state energy is always attained by a simple 
spiral with the wave vector Q that minimizes J(q), the 
Fourier transform of of the exchange parameters Jij (the 
LK theorem) . The energy per spin of a spiral with wave 
vector q is just J(q). The spiral concept was introduced 
in 1959. [1, H, This is relevant to the case discussed 
in [l|, the lattice considered being a Bravais lattice (one 
spin per primitive unit cell). The proof of the theorem 
was via the Luttinger-Tisza method 0, H, [1] ; please see 
the recent review [10| . 

The case under consideration [l| essentially involves 
spins on a square lattice, the 4 n.n. interactions being 




FIG. 1: Unit cell showing the exchange interactions between 
the Mn ions. 



equal (Ji), but, due to the internal oxygen structure, the 
n.n.n interactions J2,J3 are different. The unit cell is 
shown in Fig. 1, where the J's are indicated (the oxygens 
are not shown). Taking the x and y axes respectively, 
pointing to the right and up with respect to the figure, 
one readily finds 

J(q)/2 = Ji{co8qx+cosqy)+J2COs{qj:+qy)+J3COs{qx-qy). 

(2) 

In [ij, Ji = —\Ji\ (ferromagnetic n.n. interactions), J2 > 
and J3 is neglected, so that 

J(q)/(2| Ji|) = -{cosqr, + cosqy) + ■y cos{q^ + qy), (3) 

where 7 = J2/I Ji|. We first consider q in the (1,1) direc- 
tion, q = p{x + y), so, putting | Ji| = 1, 

J(q)/2 = -2cosp-|-7Cos(2p). (4) 

The stationary points of this are given by sinp = and 

cosp = 1/(27) = cospo; (5) 

the latter occurs only if 27 > 1 (the wave vector is real.) 
The energies at q = and at pq{x + y) = Q are, from @ 



2 





FIG. 2: upper: Energy of the ground and other states, lower: 
Component po of the spiral wave vector. 



and dH), 



J(0)/2 



-2 + 7 
1 



1 X 1 

27 



(6) 



These are equal at 7 = 1 /2 and in fact 



J(0)/2 - J(Q)/2 = 2(7 + ^ - 1) - ^(27 - 1)^ (7) 

that is, the spiral (see the LK theorem) goes below the 
ferromagnetic state as 7 increases past 1/2. This con- 
tradicts the story in [l| where the ferromagnetic state is 
indicated to be the ground state all the way to 7 — 1. 
Furthermore, the ground state is not the 'up-up-down- 
down" configuration when the ferromagnetic state is first 
destabilized, as claimed in but rather a spiral whose 
wave vector increases continuously from zero at 7 = 1/2. 
Interestingly, as 7 — > cxD,po ~^ t^/'^i which gives the spi- 
ral with a 90° turn angle; and as pointed out in [2|, this 
is degenerate with the "up-up-down-down" state. These 
results are summarized in Fig. 2. 



To complete the solution to the ground state of this 
model, we need to check that the lowest J(q) that we 
found (under the restriction q in the (1,1) direction), is 
lowest over all q in the Brillouin zone. We did this for 
two values of 7, namely, 0.4 and 0.6 (sample values in 
the ferromagnetic region and the spiral region, respec- 
tively). Indeed the (1,1) solution is the minimum over 
all q. We consider this sufficient, although it is easy to 
check further values of 7. Lest one might think that we 
have merely found the lowest spiral (the ferromagnet is 
the q = spiral), and hence ask whether there are other 
lower-energy states with more than one wave vector (and 
its negative) in their Fourier representation, we note that 
because of the place of the calculation in the framework 
of the Lutttinger-Tisza method and the LK theorem, the 
conclusion is far stronger: there is no state with a lower 
energy than the solution given. [ll| The finite-T phase 
diagram in [1;], obtained within the mean- field approx- 
imation, shows very complex behavior, numerous long- 
range-ordered phases (the "Devil's flower"). However, it 
was shown by Lyons |12| . that for a Bravais lattice with 
T-independent exchange parameters (the case considered 
in [l[), the phase boundary is vertical: a spiral at the or- 
dering temperature remains a spiral with the same wave 
vector for all T. That is, all the sweet complexity shown 
in the phase diagram of [l| is not a property of their 
assumed Hamiltonian. 



Thus the question remains, what is the source of the 
E- or "up- up-down-down" phase found in HoMnOa? It 
seems to us that the answer should be found in the stan- 
dard fundamental model of localized electrons in this in- 
sulator, as set forth by P.W. Anderson [l^, and charac- 
terized qualitatively for the Heisenberg interactions that 
are a consequence [ll] by the Goodenough-Kanamori 
rules (see[l3|)- That will include Heisenberg terms with 
more general interactions than assumed so far, as well 
as anisotropic terms and isotropic higher-spin terms, e.g. 
(Si ■ 82)^, etc. One important aspect of anisotropy is that 
it can remove the degeneracy of the 90° spiral and the 
E-state (this observation renders incorrect the argument 
in leading to their conclusion that the "classical spin 
model is not suitable to explore in a single framework the 
many phases found in RMnOa manganites" ) . 
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